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Abstract. Accretion powered pulsars exhibit a variety of lightcurves. In this paper we propose to classify the 
observed lightcurves as single or double pulsed.We analyze the lightcurves of 86 accretion powered pulsars and 
assign them to these classes. We present three datasets: first in which the classification can be easily done, second 
for which the classification is more difficult and not certain, and third for which we were unable to classify the 
pulsar because of lack of published data. We analyze a simple model in which the angles between the magnetic and 
rotation axis /3, and between the rotation axis and the line of sight 0 are random, and show that it is inconsistent 
with the data. We also present a model in which the angle between the magnetic axis and rotation axis is restricted 
and compare it with the data. This leads to an upper limit on the angle /3 < 40°. We conclude that there must 
be a mechanism that leads to alignment of the magnetic and spin axis in X-ray pulsars. 

Key words, neutron stars - X-rays 


1. The lightcurves of accretion powered X-ray 
pulsars 


Accreting neutron stars were discovered nearly more than 
30 years ago (Shklovsky 1967), with Cen X-3 the first one 
discovered which showed pulsations ( piacconi et al. 1971). 
Currently we know nearly one hundred accreting neutron 
stars, and in more than eighty of them pulsations were 
identified ( Liu et al. 200C , 2001a ). Accreting neutron stars 
in binaries exhibit a wide range of X-ray light curves. They 
vary as a function of the photon energy, and moreover in 
the transient sources the pulse shapes change with the 
variation of the luminosity. 


The pulse period in accreting sources is identified with 
the rotation of a magnetized star. As the matter from the 
companion star falls onto the neutron star it is channeled 
onto the magnetic poles by the strong magnetic field of 
the neutron star. Thus the polar caps are the places where 
most of the emission takes place. Several theoretical mod¬ 
els of radiation of magnetized accretion powered neutron 
stars have been proposed. 


The magnetized radiative transfer is solved using a dif¬ 
ference scheme (Meszaros & Nagel 1985; Bulik et al. 1992D, 


or using Monte Carlo scheme (Lamb et al 


et al. 1995), for a review see Meszaros (1992). The ra¬ 


199C|; Isenberg 


diation from each polar cap produces one pulse in the 


light curve. Depending on the emission cap physics and 
the strength of the magnetic field each pulse may have 
some additional structure. When during the rotation of 
an accreting pulsar we see two polar caps the lightcurve 
should exhibit two distinct pulses (peaks), and if only one 
cap is seen then the lightcurve is single peaked. Each pul¬ 
sar can be classified as single or double peaked depending 
on the location of the observer and the inclination be¬ 
tween the rotation and the magnetic axes. In this paper 
we classify the observed pulsars into two classes: first the 
single pulsed light curves as these for which the geometry 
is such that only one cap is seen, and the double pulsed 
light curves as those for which two caps are seen. At low 
X-ray energies this distinction may not be clear, however 
in the high energy band we clearly distinguish single or 
double pulses, and t here are no tri p le (or multiple) pulse d 
lightcurves, see e.g. Nagase ( 1989 ); Bildsten et al. ( 1997 ). 
By high energy band we mean at least above lOkeV, or 
in the pulsars where cyclotron lines were found, above the 
cyclotron lines. 


The paper is organized as follows: in section 2 we sum¬ 
marize the observations, in section 3 we analyze the ex¬ 
pected shapes and expected fraction of single and double 
peaked lightcurves, in section 4 we discuss the results. We 
summarize the conclusions in section 5. 
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2. Pulsar observations 

In this section we review the existing observations of X- 
ray pulsars and attempt to classify their lightcurves ac¬ 
cording to the number of pulses in their lightcurves. We 
summarize our classification in Tables 1, 2, and 3. For 
each pulsar we have searched the literature for X-ray ob¬ 
servations, and verified if the pulsar was observed in the 
band above lOkeV. The pulsars with such observations 
are listed in Table There are 47 entries in the Table out 
of which 21 have been classified as double and 26 as sin¬ 
gle. Additionally we have listed 25 pulsars in Table ^ for 
which no lightcurve above lOkeV was available. There are 
3 pulsars with likely double light curves and 22 with single 
in Table |^. We see that in each table there are more single 
pulsed than double pulsed objects. For completion we also 
list 14 pulsars in Table |4| for which no information on the 
type of the lightcurve was available. We define three sets 
of data corresponding to the three Tables: dataset with 
the data of Table |^, dataset S 2 with the combined data of 
Tables || and ||, and dataset 53 with all the data, and the 
entries of Table ^ are assumed to be double peaked. 


3. The expected pulsars light curves 

Let us assume that an accreting neutron star has a pure 
centered dipole field. In such case the geometry of accre¬ 
tion powered pulsars is characterized by two angles: 6 - 
the angle between the rotation axis and direction to the 
observer, and /3 the angle between the magnetic and rota¬ 
tion axes. The parameter space when only one cap is seen 
is determined by the condition: 


9 <90-f3-Ag, 


( 1 ) 


where Ag is the gravitational light deflection of a ray 
emitted parallel to the surface of the neutron star. Here 
we assume that the magnetic beaming is not very strong 
and that there are X-ray photons emitted parallel to the 
surface of the neutron star. Ignoring the gravitational 
light bending effects, i.e. in the classical Newtonian case, 
Ag = 0°, the parameter space for seeing just one cap is 
shown as the region below the solid line in Figure]^. Thus, 
if the angles /3 and 9 are chosen randomly the fraction of 
single pulsed objects should be 1 — 7r/4 ss 0.214. The re¬ 
gion for the single peaked light curves would be smaller 
if the radiation were strongly beamed. However, in this 
case the double peaked light curves would exhibit narrow 
pulses, contrar y to w hat is observed - see e.g. Figure 7 in 

(S- 

For a set of data, like pulsars, with a well defined prop¬ 
erty, i.e. either single or double pulsed, and a given proba¬ 
bility p that each pulsar will belong to one or another class 
the probability of observing n single peaked lightcurves 
and m double peaked in an observation of randomly cho¬ 
sen n -|- m objects is given by the binomial distribution: 


Bildsten et al. 


P — 

n,.m. — 


(n -I- m)\ 


p^a-pY 


( 2 ) 


and the probability of obtaining M or less double peaked 
objects in a sample containing N source is a sum 


M 


'Pnm = ^ Pn-t 


m—1 


( 3 ) 


We can test the consistency of the datasets defined above 
with the simplest model: random inclination of the mag¬ 
netic axis to the rotation axis /3, random direction of the 
line of sight 9, and no gravitational light deflection. For 
the dataset Si we obtain the probability 2.66 x 10“^ that 
there are 21 or less double pulsars in a set of 47 observed 
objects. Applying this test to the dataset ^2 we obtained 
the probability of only 9.33 x 10“^^, and for the dataset 
S 3 this probability is 1.76 x 10“^^. Thus we see that the 
data are inconsistent with this simple model. This is due 
to the fact that the parameter space for viewing single 
pulsed objects is smaller than the double pulsed objects. 
Including the effects of the gravitational light deflection 
only reduces the probability of seeing one pulsed objects 
and makes the values of the above probabilities drop even 
further. We present the probability of seeing one polar cap 
as a function of Ag in the left panel of Figure ^ Another 
possibility to consider is that the pulsar beam is not wide 
but strongly collimated, and the beam has an opening an¬ 
gle of a. In this case for a given (3 there is only a small 
range of angles (3 — a<9<(3-\-aiov which the pulsar is 
visible at all. The restriction on the parameter space for 
small a = 10° is shown in Figure as the shaded region. 
The probability of seeing only one accretion cap as a func¬ 
tion of the angle a is shown in the right panel of Figure |[ 
Another possibility to consider is to allow for non aligned 
dipole, i.e. to assume that the two polar caps are not an¬ 
tipodal but can be placed anywhere on the stellar surface. 
A straightforward calculation neglecting the gravitational 
light bending effects shows that the probability of seeing 
just one polar cap drops in this case to 17%. Thus the 
probabilities computed above in the Newtonian case are 
conservative and inclusion of the gravitational light bend¬ 
ing and beaming make the result even more significant. 

The incompatibility of the simplest model with the 
data prompts us to consider a more general model. We 
assume first that the distribution of the inclinations of 
the magnetic axes to the rotation axis is restricted, i.e. 
the angles /3 are limited to the range between 0 and /3max- 
This corresponds to restricting the parameter space to the 
region above a horizontal line at cos/3niax in Figure 0. For 
each value of /3max we can compute the probability of ob¬ 
serving the given set of data using equation (i- We plot 
this probability as a function of /3max in Figure The 
standard model corresponds to the case /3max = 90° and 
the probability of seeing the data rises to nearly unity for 
the angles /3max < 40° for the datasets 5i and S 3 , and 
/3max < 30° for the dataset S 2 - The probabilities of see¬ 
ing one cap increases with decreasing /^max, and therefore 
for its intermediate values the probability of seeing the 
data is rising to nearly unity. For the very small values of 
/3max one expects to see nearly all pulsars as single pulsed. 
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Table 1. X-ray pulsars in group one, i.e. these that are easily classifiable as single an double pulsed. Part 1 


System 

P . 

J- spin 

Type 

Comments 

SAX J1808.4-3658 

0.002494 

Single 


EXO 1024.0-5732 

0.061 

Double 


GRO J1744-28 

0.467 

Single 

The Bursting Pulsar 

SMC X-1 

0.717 

Double 


Her X-1 

1.24 

Double 

2nd pulse weak 

SMC X-2 

2.37 

Single 


4U 0115-163 

3.61 

Double 


V0332+53 

4.37 

Single 


GRO J1750-27 

4.45 

Single 


Cen X-3 

4.82 

Double 

2nd pulse weak 

IE 2259-1586 

6.98 

Double 


4U 1626-67 

7.66 

Single 


GS 0834-430 

12.3 

Double 


LMC X-4 

13.4 

Single 


J00521-7319 

15.3 

Double 


XTE J1946+274 

15.8 

Double 


2S 1417-624 

17.6 

Double 


GRO J1948+32 

18.7 

Single 


RX J0117.6-7330 

22.0 

Double 


GS 1843+00 

29.5 

Double 


J0112.2-7317 

31.03 

Single 


RX J0812.4-3114 

31.9 

Single 

accretion flow eclipse 

OAO 1657-415 

37.7 

Single 


EXO 2030+375 

41.7 

Double 

broad peak with two maxima 

ISAX J0054.9-7226 

59 

Single 


Cep X-4 

66.2 

Double 


XTE J1906+09 

89.2 

Single 

double pk below 20keV 

GRO J1008-57 

93.5 

Single 


GS 1843-024 

94.8 

Single 


4U 0728-25 

103.2 

Double 


A0535+26 

105 

Double 

accretion flow eclipse 

Set X-1 

111 

Double 

broad peak two maxima 

GX 1+4 

120 

Single 

accretion flow eclipse 

GRO J2058+42 

198 

Single 


Vela X-1 

283 

Double 
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Table 2. X-ray pulsars in group one, i.e. these that are easily classifiable as single an double pulsed. Part 2 


System Pspin Type Comments 

4U 1145-619 292 Sin g le 


IE 1145.1-6141 

297 

J0103.2-7209 

345.2 

SAX J2103.5-b4545 

358.6 

XTE J1855-026 

361.1 

A 1118-616 

405 

GPS 1722-363 

414 

4U 1907-b09 

438 

4U 1538-52 

530 

GX 301-2 

681 

X Per 

835 

RX J0146.9-t6121 

1413 


Single 

single 

Single 

Single 

Single 

Single 

Double 2nd pulse weak 

Double 

Double 

Single 

Single 


References 


^ jaiasw at el al.^ (1997); Bildsten et al, 
( [1997D 

kdievenev et al.* (1992) 

Israel et al.^ (200C) 

*Hulleiiiaji el at? (1998) 
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' Ooe et al. l ( |l994|) 
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tjild&ten et al.^ (1997) 
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* Rubba &. Waiwick* (1989); Robba et al 
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*Mei ' eglietti et al -* (2000) 




Fig. 1. Parameter space for viewing accreting X-ray pul¬ 
sars. The region below the solid line corresponds to vis¬ 
ibility of two polar caps with gravitational light bending 
neglected. The region below the dashed line corresponds 
to visibility of two caps with the gravitational light de¬ 
flection Ag = 30°. The shaded region is the part of the 
parameter space in which pulses would be detectable if 
radiation is strongly beamed. The horizontal dotted line 
corresponds to the angle between the rotation and mag¬ 
netic axes f3 = 40°. 

Since each dataset contains some double pulsed objects 
the probability of seeing the data drops when /3max ^ 0. 


Fig. 2. Left panel: the probability of seeing only one polar 
cap as a function of the maximal gravitational light de¬ 
flection. The non-relativistic case corresponds to Ag = 0. 
Right panel; the probability of seeing only one cap as a 
function of the beaming angle a. 

A more general model taking into account the general 
relativistic effects includes two parameters: the maximal 
angle between the magnetic and rotation axis /9max and 
the maximum gravitational light deflection angle Ag. For 
such a model we compute the probability of seeing a pul¬ 
sar with only one polar cap as a function of these two 
angles. In Figure ^ we present the regions containing the 
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Table 3. X-ray pulsars in group two, i.e. with no pulse shape above lOkeV or were difficult to classify. 


System 

P . 

■L spin 

Type 

Comments 

SAX J0635+0533 

0.0338 

Single 

below 10 keV 

RX J0059.2-7138 

2.76 

Single 

ROSAT obs 

RX J0502.9-6626 

4.06 

Single 

ROSAT obs 

RX J1838.4-0301 

5.45 

Single 

ROSAT obs 

IE 1048-593 

6.44 

Single 

Einstein obs 

AX J1845.0-0300 

7 

Single 

ASCA obs 

4U 0142-b614 

8.69 

Double 

obs below 10 keV 

2E 0050.1-7247 

8.88 

Single 

ROSAT obs 

2S 1553-54 

9.26 

Double 

SAS-3 obs below llkeV 

AX J0049-732 

9.31 

Single 

ASCA obs 

IRXS 170849.0-400910 

11.0 

Single 

ASCA obs; could be an AnXP 

RX J0648.1-4419 

13.2 

Single 

ROSAT obs 

RX J0529.8-6556 

69.5 

Single 

ROSAT obs 

AX J0049-729 

74.8 

Single 

ASCA obs 

SAX J0544.1-710 

96.1 

Single 

obs below lOkeV 

ISAX J1324.4-6200 

170 

Single 

SAX obs below 20keV 

RX J0440.9-t4431 

202.5 

Single 

obs below 20 keV 

AX J1749.2-2725 

220 

Single 

ASCA obs 

GX 304-1 

272 

Single 

SAS-3 obs 

AX J0058-7203 

280 

Double 

ASCA obs 

AX J0051-733 

323 

Single 

ASCA obs 

ISAX J1452.8-5949 

437 

Single 

obs MECS below lOkeV 

AX J170006-4157 

714.5 

Single 

ASCA obs 

RX J1037.5-564 

860 

Single 

pulse shape below 20 keV 

2S 0114+650 

10008 

Single 

ASM obs below 12 keV 
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models that are consistent with the data with the proba¬ 
bility higher than 40%, 10%, and 1%. The non relativistic 
results correspond to the line Ag = 0 in Figure Since 
the probability of seeing only one polar cap decreases with 
increasing Ag, we also set a limit on its maximum value. 
The region corresponding to the 1% agreement limit for 
the dataset can be expressed as Ag + 0.63/3i„ax < 33°. 
Thus if there were no gravitational light deflection ef¬ 
fects we find that the inclination angle of the magnetic 
axis to the rotation axis should be smaller than ss 50°. 
For the case with a moderate gravitational light bending 
Ag = 20° the magnetic axis inclination must be smaller 
than /3max < 20°. We also see that the amount of gravita¬ 
tional light bending cannot exceed 33° at the 99% confi¬ 
dence level. The limits are even stronger when we take into 


account dataset S 2 , however for the set S 3 the constraints 
are similar to these obtained for the dataset 5i. 

4. Discussion 

To discuss the limit on Ag obtained above we will consider 
two cases: (i) emission from a polar cap, and (ii) emission 
from the accretion column. In the first case the upper limit 
on the Ag obtained above can be interpreted as a limit on 
the compactness of the star. We plot the amount of gravi¬ 
tational light deflection as a function of the emission angle 
for a few values of the ratio of the neutron star radius to 
the gravitational radius r/vg in Figure The maximal 
gravitational light deflection of about 30° corresponds to 
the case with ratio r/vg 3 for the photons emitted par¬ 
allel to the surface. The limit on Ag can be understood as 
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Table 4. X-ray pulsar in group three, for which no reliable pulse shapes were found. 


System Papin Type Comments 


A0538-67 

0.069 

7 

no Ic published 

J0501.1-7211 

3.343 

? 

no Ic published 

RX J0720.4-3125 

8.38 

7 

isolated ns accreting from ISM? 

EXO 053109-6609.2 

13.7 

7 

EXOSAT detection 

J0053.8-7226 

46.6 

7 

no ic published 

X 1839-04 

81.1 

7 

no Ic published 

J0051-722 

91.1 

7 


AX 1820.5-1434 

152 

7 

weak pulsations 

J0054-720 

169.3 

7 

no Ic published 

3U 1239-599 

191 

7 

Ariel obs, no Ic published 

XTE J1858+034 

221 

7 

no Ic published 

4U 2206 -k 54 

392 

7 

no Ic published 

IWGA J1958.2-k3232 

730 

7 

only ROSAT Ic 

AX J0049.5-7323 

755 

7 

weak pulsations 
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Fig. 4. Each panel presents the region containing the models that are consistent with the data at the level of 60%, 
90%, and 99%. The three panels correspond to the three datasets. 


the upper limit on the compactness of the neutron stars. 
In Figure we plot the ratio r/vg as a function of the 
stellar gravitational mass for a few representative equa¬ 
tions of state. The ratio r/vg decreases with increasing 
mass. The upper limit on the ratio r/rg corresponds for 
a given equation of state to the lower limit on the mass 
of the accreting neutron stars. Masses of neutron stars in 
some accreting objects have been estimated for from the 
observations of the binary motion, like e.g. in the case of 
4U1538-52 (Pakull et al. 1983; Reynolds et al. 1992| ) and 
they all seem to be consistent with the fiducial 1.4 M 0 
with some uncertainty. The general consensus seems to be 
that the masses of all neutron stars are consistent with the 
mean value of w 1.4M^ , and this is also the value foun d 
for the binary pulsars (Thorsett fc Chakrabarty 1999|), 


however in LMXBs this range may be much wider and 
reach above 2 Mq Bulik et al. ( 2000 ). If the neutron stars 
have masses around 1.4 Mq and the beams are wide then 
our results point to the very stiff equations of state like 
e.g. SDB (Sahu et al. 1993). 


A second possibility is that we see the emission from 
the accretion column with the typical properties of the fan 
beam and that this emission is responsible for the observed 
flux at angles with respect to the normal above 90 degrees. 
Since we limit the width of the beam to Ag < 33°. In 
this case our limit for the non-relativistic case means that 
the emission from the column extends no higher than to 
« 0.2r. If the effects of the gravitational light bending 
around a neutron star are included the maximum height 
allowed is even smaller . Thus our results argue for the 
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Emission angle [deg] 


Fig. 3. The probability of observing the data in models 
parameterized by the the maximum angle between the 
magnetic and rotation axis f3max- The solid line corre¬ 
sponds to the data set 5i, and the short-dashed line cor¬ 
responds to the data set <S 2 . The long-dashed line corre¬ 
sponds to the data set 53 . 


Fig. 5. The value of the gravitational light bending as 
a function of the emission angle, measured from the ra¬ 
dial direction for photons emitted at different initial radii. 
The maximum gravitational light deflection is the value 
of deflection for the initial angle of 90°, i.e. parallel to the 
surface. 


polar emission and against signihcant emission from the 
accretion column. 

We have obtained a a limit on the possible inclination 
angles of the magnetic axis with respect to the rotation 
axis f3max < 50°. A hint of such distribution was already 


shown by Leahy ( 1991 ). He analyzed the lightcurves of 20 
pulsars using a simple model with the beam dependence 
on the emission angle as cos^ (f) and cos'* </> where 4 > is the 
direction to the normal of the surface. For each pulsar he 
found the angles 9 and [3 however he fits were degenerate 
in these angles, so that it was not clear a priori which an¬ 
gle is which. The resulting distribution was inconsistent 
with the assumption that both angles are random, and 
Leahy (|1991) concluded that the mean angle between the 


magnetic and rotation axis is as small 15°. Bulik et al. 


( 1995| ) modeled lightcurves of 4U 1538-52 and Vela X-1 


observed by Ginga, and obtained the inclination angles 


below 50° for their hts with best significance. Kraus et al. 


(1996) analyzed the pulse profiles of Cen X-3 and recon¬ 
structed the geometry and beam pattern for this object. 



Fig. 6. The ratio of the stellar radius to the gravitational 
radius as a function of mass for several neutron star equa¬ 
tions of state. The labels are BPAL12 - Prakash et al 


(19971), SBD - Sahu et al. (1993), BBl and BB2 - Baldc 


They found that the most likely inclination angles of the al.| ( |I^ , while SSI and SS2 corr espond to the MIT 
magneticaxis^ axis are 20°. Bag model of strange quark matter (Witten 1984) with 

Cemeljic & Bulik ( 1998 ) have noticed that the light different densities at zero pressure, 

curves of some accretion powered pulsars exhibit sud¬ 
den dips. Such dips can be interpreted as accretion flow 
eclipses, i.e. be caused by the absorption of the X-ray flux 
when the accretion column passes the line of sight. Such 
phenomena have been noted for at least three objects: 


A0535-262 (Cemeljic & Bulik 1998), GXl-l-4 (Giles et al. 
2 OOOI), and LS992 (|Reig fc Roche||199'^. 


If such eclipses happen when the accretion column 
passes through the line of sight then they should take place 
when \9 — (3\ < 7 , where 7 is the angular width of the ac¬ 
cretion column. The region corresponding to 7 = 10° is 
shaded in Figure |^. We see that such eclipses should fre- 
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quently take place in double pulsed light curves, if the 
magnetic axis is randomly oriented with respect to the ro¬ 
tation axis. Moreover within the single pulsed light curves 
a signihcant fraction, - approximately 20-30% should ex¬ 
hibit such dips, and this fraction increases as the amount 
of gravitational light bending increases. 

Yet the only three cases when such dips have been 
found are in single pulsed light curves. The lightcurve of 
A0535-262 seems to be double peaked, however Cemeljic 
& Buli^ ( |1998[ ) pointed out that it can be interpreted as 
a single pulse divided by the accretion flow eclipse. This 
lack of double peaked pulsars with accretion flow eclipses 
(dips) may mean that the region of parameter space for 
double peaked light curves with dips is either empty or 
scarcely populated! Again we are drawn to the conclusion 
that the magnetic axis needs to be aligned with the rota¬ 
tion axis. If the magnetic axis inclination is restricted to 
the region (3max < 50°, than one expects that the accre¬ 
tion flow eclipses to be observed only in the one pulsed 
lightcurves, as is the case. While this by itself is not a sig¬ 
nificant result, the identification of accretion flow eclipses 
is consistent with the conclusion that the range of allowed 
values of (3 is limited. 

One should mention here that there may possibly be 
a selection bias against seeing double pulsed objects. In 
the case of nearly perpendicular rotators the two pulses 
will be similar. It is therefore possible that the lack of 
double pulsed objects is due to misidentiflcation of the 
pulse period which should be two times longer as reported. 
The fraction of such misidentified objects depends on the 
exact shape of the beam. For wide uncollimated beams the 
two pulses look similar for a range of angles f3 as large as 
the width of the beams. However, our limit is (3max < 50°, 
so the beams would have to be very wide. Additionally, the 
beam shape would have to be very fine tuned in order to 
compensate the viewing geometry effects over such large 
range of angles. 

The sample of accreting pulsars presented here is dom¬ 
inated by transient Be/X-ray sytems. One should note 
that for the disk fed persistent sources (SMC X-1, SMC 
X-2, Cen X-3, LMC X-4, RX J0648.1-4419) there are two 
double peaked objects out of five sources, still below the 
expected ~ 79%, but this is not a significant discrep¬ 
ancy. In the case of wind fed persistent sources (Vela X- 
1, IE 1145.1-6141, 4U1907-k09, 4U1583-53, GX 301-2, 2S 
0114-1-650) there are four double peaked out of six objects. 
However, because of the small number of these sources it 
is impossible to obtain any significant conclusion just from 
statistical analysis. However, as already mentioned above 
detailed modeling of two wind fed pulsars (4U 1538-52 
and Vela X-1) and one disk fed (Cen X-3) resulted in low 
values of (3. 


5. Conclusions 

We analyzed the shapes of accretion powered pulsar light 
curves, looking for the expected fraction of single and dou¬ 
ble peaked light curves, as well as these exhibiting accre¬ 


tion flow eclipses. We find that the expected number of sin¬ 
gle peaked lightcurves is small; certainly less than about 
one in five lightcurves should be single peaked. On the 
other hand the observations show that quite a large frac¬ 
tion of accreting pulsars have single peaked lightcurves 
in hard X-rays. In order to explain that we are drawn 
to conclude that the magnetic axes in accretion pow¬ 
ered pulsars tend to be aligned with the rotation axis. 
Such alignment decreases the expected number of double 
peaked lightcurves, as seen in Figure]^. We have compared 
the models parametrized by the maximal inclination of 
the magnetic to the rotation axis Pmax and the value of 
the maximum gravitational light deflection Ag with the 
data. We obtain upper limits at the 99% confidence level 
/3max < 50° and Ag < 33°. 

Our results are based on the classification of the neu¬ 
tron star lightcurves as single or double pulsed. In order 
to avoid possible misclassification we have attempted to 
use the data above 10 keV, or above the cyclotron line 
wherever possible. Moreover, we have attempted to clas¬ 
sify as many pulsars as double when we were in doubt. 
In order to rid of the possible systematic effects we have 
repeated the analysis for the three datasets, and in the 
third dataset we included all the unclassiflable pulsars as 
the double peaked objects to be conservative. The deficit 
of double pulsed objects is significant for every dataset 
that we analyzed. Moreover we note that the paucity of 
the double peaked pulsars is so strong that our results still 
hold even if there are still a few misidentifications in our 
analysis. We have also shown that if the pulsar beam is 
strongly collimated the deficit of double peaked objects 
becomes more significant. 
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